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Abstract. In this article, we study the numerical approximation of stochastic differ- 
ential equations driven by a multidimensional fractional Brownian motion (ffim) with 
Hurst parameter greater than 1/3. We introduce an implcmentable scheme for these 
equations, which is based on a second order Taylor expansion, where the usual Levy 
area terms are replaced by products of increments of the driving fBm. The convergence 
of our scheme is shown by means of a combination of rough paths techniques and error 
bounds for the discretisation of the Levy area terms. 



1. Introduction and Main Results 

Fractional Brownian motion (fBm in short for the remainder of the article) is a nat- 
ural generalisation of the usual Brownian motion, insofar as it is defined as a centered 
Gaussian process B = {B t ; t G M+} with continuous sample paths, whose increments 
(SB) st := B t — B s , s,t G M+ are characterised by their variance E[(SB)l t ] = \t — s\ 2H . 
Here the parameter H G (0, 1), which is called Hurst parameter, governs in particular the 
Holder regularity of the sample paths of B by a standard application of Kolmogorov's 
criterion: fBm has Holder continuous sample paths of order A for all A < H. The par- 
ticular case H = 1/2 corresponds to the usual Brownian motion, so the cases H ^ 1/2 
are a natural extension of the classical situation, allowing e.g. any prescribed Holder 
regularity of the driving process. Moreover, fBm is if -self similar, i.e. for any c > the 
process {c H B t / c ; t G M + } is again a fBm, and also has stationarity increments, that is 
for any h > the process {B t+ h — B h ; t G 1R + } is a fBm. 

These properties (partially) explain why stochastic equations driven by fBm have 
received considerable attention during the last two decades. Indeed, many physical sys- 
tems seem to be governed by a Gaussian noise with different properties than classical 
Brownian motion. Fractional Brownian motion as driving noise is used e.g. in electrical 
engineering [121 [13], or biophysics [5j|23j[3l]. Moreover, after some controversial discus- 
sions (see [3] for a summary of the early developments) fBm has established itself also in 
financial modelling, see e.g. [TT], [2] . For empirical studies of fractional Brownian motion 
in finance see e.g. [8j |39j [7]. All these situations lead to different kind of stochastic 
differential equations (SDEs), whose simplest prototype can be formally written as 



Y 



/ v®(Y u )dB%\ tG[0,T], aeR d , (1) 
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where a = (a^\ . . . ,a^) is a smooth enough function from M. d to M dxm and B = 
(B^\ . . . , B^) is a m-dimensional fBm with Hurst parameter H > 1/3. 

At a mathematical level, fractional differential equations of type (CEJ) are typically 
handled (for H ^ 1/2) by pathwise or semi-pathwise methods. Indeed for H > 1/2, the 
integrals J* (Y u ) dBu\ i = l,...,m, in ([1]) can be defined using Young integration 
or fractional calculus tools, and these methods also yield the existence of a unique 
solution, see e.g. [331 SO]- When 1/4 < H < 1/2, the existence and uniqueness result 
for equation ([1]) can be seen as the canonical example of an application of the rough 
paths theory. The reader is referred to [16], [25] for the original version of the rough paths 
theory, and to [TS] for a (slightly) simpler algebraic setting which will be used in the 
current article. In the particular case 1/3 < H < 1/2, the rough path machinery can 
be summarised very briefly as follows: assume that our driving signal B allows to define 
iterated integrals with respect to itself. Then one can define and solve equation ([T]) in a 
reasonable class of processes. 

Once SDEs driven by fBm are solved, it is quite natural (as in the case of SDEs 
driven by the usual Brownian motion) to study the stochastic processes they define. 
However, even if some progress has been made in this direction, e.g. concerning the law 
of the solution [U [H [25] or its ergodic properties [2U] , the picture here is far from being 
complete. Moreover, explicit solutions of stochastic differential equations driven by fBm 
are rarely known, as in the case of SDEs driven by classical Brownian motion. Thus one 
has to rely on numerical methods for the simulation of these equations. 

So far, some numerical schemes for equations like ([[]) have already been studied in the 
literature. In the following, we consider uniform grids of the form {tk = kT/n; < k < 
n} for a fixed T > 0. The simplest approximation method is the Euler scheme defined 
by 

Y n = a, 

m 

y t : +1 = ^:+E * w (i2)*sSL> * = o, . . . , n - 1. 

i=l 

For H > 1/2, the Euler scheme converges to the solution of the SDE ([!)). See e.g. 
in [25], where an almost sure convergence rate n -( 2H ~ 1 )+ £ with e > arbitrarily small is 
established. A detailed analysis of the one-dimensional case is given in [28], where the 
exact convergence rate n~ 2H+1 and the asymptotic error distribution are derived. 

However, the Euler scheme is not appropriate to approximate SDEs driven by fBm 
when 1/3 < H < 1/2. This is easily illustrated by the following one-dimensional exam- 
ple, in which B denotes a one-dimensional fBm: consider the equation 

dY t = Y t dB u te[0,l], Y = l, 

whose exact solution is 

Y t = exp(B t ), t e [0,1]. 

The Euler approximation for this equation at the final time point t — 1 can be written 

as 

n-l 

Yr = II^ 1 + (SB)k/n,(k + l)/n). 
k=0 
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So for neN sufficiently large and using a Taylor expansion, we have 

n—1 ^ n— 1 

*T = eX P (E 1 ^ 1 + ( SB )k/n,(k+l)/n)j = exp (#i - - ^ |(<JS)fc/ n ,(fc+l)/n| 2 + Pn) , 
fc=0 A;=0 

where lim^oo p n =' for H > 1/3. Now it is well known that 



n-l 

^ |(^-B)fe/n,(fe+l)/n| 2 
fc=0 



OG 



for if < 1/2 as n — > oo, which implies that lim^oo Y™ =' 0. This is obviously incompat- 
ible with a convergence towards Y\ = exp (Si). In the case H = 1/2 this phenomenon 
is also well known: here the Euler scheme converges to the Ito solution and not to the 
Stratonovich solution of SDE (pQ). 

To obtain a convergent numerical method Davie proposed in [H] a scheme of Milstein 
type. For this, assume that all iterated integrals of B with respect to itself are collected 
into a m x m matrix B 2 , i.e. set 



B 2 st (t,j)= dB®dB®, 0<s<t<T, l<i,j<m. 

J s J s 

The matrix B 2 (respectively its elements) is (are) usually called Levy area. Davie's 
scheme is then given by 

Yo = a, (2) 

m m 

K + r = K + J2^( Y SB ti +1 + £ Kt k Jhj), k = o, . . . ,n - 1, 

i=l i,j = l 

with the differential operator X)W = J^ =1 a\ ^ d Xl . (Recall that we use the notation 
5Bgt = Bjr — Bs for s,t G [0, T].) This scheme is shown to be convergent as long as 
H > 1/3 in j9], with an almost sure convergence rate of n -( 3H ~ 1 )+ e for e > arbitrarily 
small. This result has then been extended in [16] to an abstract rough path with arbitrary 
regularity, under further assumptions on the higher order iterated integral of the driving 
signal. 

As the classical Milstein scheme for SDEs driven by Brownian motion, the Milstein- 
type scheme fl2]) is in general not a directly implementable method. Indeed, unless the 
commutativity condition 

V®<t® = V®<t®, i,j = l,...,m, 

holds, the simulation of the iterated integrals B 2 fetfe+i (i, j) is necessary. However, the law 
of these integrals is unknown, so that they can not be simulated directly and have to be 
approximated. 

In this article we replace the iterated integrals by a simple product of increments, i.e. 
we use the approximation 

B 2 (i,j)n-6Bi% SB® . (3) 

t k t k+ i\ > J/ r> ifcifc+i tfctfc+i V / 
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This leads to the following simpler Milstein-type scheme: Set Z™ = a and 

m _ m 

for fc = 0, . . . , n — 1. Moreover, for t G (tfc, tfc+i), define 

Z" = Z" +^- 7 ^(5Z n ) t + , (5) 
k Tin v 7 ***fe+i v y 

i.e. if £ G [0, T] is not a discretisation point, then Z" is defined by piecewise linear 
interpolation. This scheme is now directly implementable and is still convergent. 



Theorem 1.1. Assume that a G C 3 (lR d ; M dxm ) is bounded with bounded derivatives. 
Let Y be the solution to equation (QJ) and Z n the Milstein approximation given by OL) 
and Moreover, let 1/3 < 7 < i/. T/ien, £/iere exists a finite and non-negative 
random variable tjh^^t such that 

\\Y - Z n \\^ T < VH ^ T • v/logR) • (6) 

for n > 1 . 



Here || • || Kj0 o,r denotes the K-Holder norm of a function / : [0, T] — > W, i.e. 

\\f\U,oo,T= snp \f{t)\+ sup lf ^~ f ] i K S)l . (7) 
te[o,T] s,te[o,T] \ l — s \ 

Remark 1.2. Note that the almost sure estimate ([6]) cannot be turned into an /^-estimate 
for || Y — Z n \\ 1)00 ,T- This is a common consequence of the use of the rough paths method, 
which exhibits non-integrable (random) constants, as a careful examination of the proof 
of Theorem 12.61 would show. See also [16] for further details. 



Our strategy to prove the above Theorem consists of two steps. First we determine 
the error between Y and its Wong-Zakai approximation 

rn „ t 

z n t =*+Y. ail) ^ dB u ,n > tE [°» T i > ae Rd > ( 8 ) 



where 



i.e. B in equation flTJ is replaced with its piecewise linear interpolation. (For a survey 
on Wong-Zakai approximations for standard SDEs see e.g. [36J.) Here, we denote the 
Levy area corresponding to B n by B n . Using the Lipschitzness of the It 6 map of Y, 
i.e. the solution of equation ([Q) depends continuously in appropriate Holder norms on 
B and the Levy-area B, and error bounds for the difference between B and B n resp. B 
and B n , we obtain 



where T is a finite and non-negative random variable. 

In the second step we analyse the difference between ~Z a and Z n . The second order 
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Taylor scheme with stepsize T/n for classical ordinary differential equations applied to 
the Wong-Zakai approximation (JSJ) gives our simplified Milstein scheme (j3]). So to obtain 
the error bound 

\\Z n - ^1 7) oo.t < rfS^T ■ Vlog(n) • n-^), 

we can proceed in a similar way as for the numerical analysis of classical ordinary 
differential equations. We first determine the one-step error and then control the error 
propagation using a global stability result with respect to the initial value for differential 
equations driven by rough paths. The latter can be considered as a substitute for 
Gronwall's lemma in this context. 

Combining both error bounds then gives Theorem 11.11 

Remark 1.3. For H = 1/2 the scheme (j2J) corresponds to the classical Milstein scheme 
for Stratonovich SDEs driven by Brownian motion, while our scheme (J4]) corresponds to 
the so called simplified Milstein scheme. See e.g. [22] . 

Remark 1.4. At the price of further computations, which are simpler than the ones in 
this article, our convergence result can be extended to an equation with drift, i.e. to 

nt m ft 

Y t = a+ / b(Y u ) du + J2 dBf, t £ [0, T] , a £ R d , 

Jo i=1 Jo 

where b : M. d — > IR d is a function and where the other coefficients satisfy the as- 
sumptions of Theorem 11.11 Indeed, the equation above can be treated like our original 
system (QQ) by adding a component B^ = t to the fractional Brownian motion. The 
additional iterated integrals of with respect to B^) for j = l,...,m are easier 
to handle than B 2 (z, j) for i,j £ {1, . . . , m}, since they are classical Riemann-Stieltjes 
integrals. For sake of conciseness we do not include the corresponding details. 

Remark 1.5. Theorem 11.11 requires a to be bounded. However, if a £ C 3 (lR ,i ; M dxm ) is 
neither bounded nor has bounded derivatives but equation (JT|) has still a unique pathwise 
solution in the sense of Theorem 12.61 below, then the assertion of Theorem 11.11 is still 
valid. This follows from a standard localisation procedure, see e.g. [21], and applies in 
particular to affine-linear coefficients. 

Remark 1.6. The error bound of Theorem 11.11 is sharp. To see this, consider the most 
simple equation 

dY t {1) = dBf\ t £ [0, T] , Y = aeR, 

for which our approximation obviously reduces to Z n = B n . Then, due to results of 
Hiisler, Piterbarg and Seleznjev ([H]) for the deviation of a Gaussian process from its 
linear approximation, one can prove that 

lim P (£(n) ■ \\Y - Z n || 7i00jT < oo) = 0, 

n— ¥oo 

if 

liminf l(n) ■ ^log(n) • n~^ H ~^ = oo. (9) 

n— >od 

For further details see Section [TBI 
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Remark 1.7. If the Wong-Zakai approximation is discretised with an arbitrary numerical 
scheme for ODEs of at least second order (e.g. Heun, Runge-Kutta 4), then the arising 
scheme for equation ([I]) satisfies the same error bound as the proposed modified Milstein 
scheme. So, the strategy of our proof is in fact an instruction for the construction of 
arbitrary implementable and convergent numerical schemes for SDEs driven by fBm. 

Remark 1.8. Instead of replacing the Levy terms in Davie's scheme by the "rough" 
approximation one could discretise these terms very finely using the results contained 
in [21], where (exact) convergence rates for approximations of the Levy area are derived. 
However, it is well known that already for SDEs driven by Brownian motion such a 
scheme is rarely efficient, if the convergence rate of the scheme is measured in terms of 
its computational cost. For a survey on the complexity of the approximation of SDEs 
driven by Brownian motion, see e.g. [27]. 

The 7- Holder norm, which appears in Theorem 11.11 since the Ito-map of Y is only 
Lipschitz in appropriate Holder norms with 1/3 < 7 < H and thus is natural in the 
rough path setting, is not typical for measuring the error of approximations to stochastic 
differential equations. A more standard criterion would be the error with respect to the 
supremum norm, i.e. 

\\Y-Z n \\^ T = sup \Y t -Z?\. 

te[o,T] 

The error (in the supremum norm) of the piecewise linear interpolation of fractional 
Brownian motion is of order ^\og{n)n~ H , see [T4] . Moreover, for the iterated inte- 
gral / T J " dB^ dBu^ the proposed Milstein-type scheme leads to the trapezoidal type 
approximation 

1 n— 1 

iE(4 1, +s£,)(4 2 i 1 -<)- 

k=0 

The L p -error for this approximation is of order n~ 2H+1 ^ 2 , see [31] . 

Based on these two findings, our guess for the rate of convergence in supremum norm 
is that 

\\Y - Z n \\^ T < VH ^ ■ V^n) ■ {n~ H + n- 2H+l ' 2 ) 
holds under the assumptions of Theorem 11.11 This conjecture is also supported by the 
numerical examples we give in Section 4. 

The remainder of this article is structured as follows: In Section [2] we recall some 
basic facts on algebraic integration and rough differential equations. The proofs of 
Theorem 11.11 and Remark 11.61 are given in Section 3 and 4. Finally, Section 5 contains 
the mentioned numerical examples. 

2. Algebraic integration and differential equations 

In this section, we recall the main concepts of algebraic integration, which will be 
essential to define the generalized integrals in our setting. Namely, we state the definition 
of the spaces of increments, of the operator 5, and its inverse called A (or sewing map 
according to the terminology of [15]). We also recall some elementary but useful algebraic 
relations on the spaces of increments. The interested reader is sent to [IB] for a complete 
account on the topic, or to [TTT [T§] for a more detailed summary. 
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2.1. Increments. The extended integral we deal with is based on the notion of incre- 
ments, together with an elementary operator 5 acting on them. 

The notion of increment can be introduced in the following way: for two arbitrary real 
numbers £2 > £\ > 0, a vector space V, and an integer k > 1, we denote by C k {\£\, £2}', V) 
the set of continuous functions g : [£i,£2] k — > V such that gt v -t k = whenever U = U+i 
for some i G {0, . . . , k — 1}. Such a function will be called a {k — 1) -increment, and 
we will set C*([£i, £2]; ^0 = U k >xC k ([£i, £2]', V). To simplify the notation, we will write 
Ck(V), if there is no ambiguity about [li,^]- 

The operator 5 is an operator acting on k- increments, and is defined as follows on 
C„{V): 

fc+i 

5 : C k (V) C k+1 (V), (Sg) tl .., k+1 = ^(-1)^..,;..,^, (10) 

i=l 

where t« means that this particular argument is omitted. Then a fundamental property 
of 5, which is easily verified, is that 55 = 0, where 55 is considered as an operator from 
C k (V) to C k+2 {V). We will denote ZC k {V) = C k {V) n Ker<5 and BC k {V) = C k (V) n lm5. 

Some simple examples of actions of 5, which will be the ones we will really use through- 
out the article, are obtained by letting g G Ci(V) and h G C 2 {V). Then, for any 

t,u,s G [£i, £2], we have 

(5g) st = g t -g s and {5h) sut = h st - h su - h ut . (11) 

Our future discussions will mainly rely on ^-increments with k = 2 or k = 3, for 
which we will use some analytical assumptions. Namely, we measure the size of these 
increments by Holder norms defined in the following way: for / G ^(V) let 

11/11, = sup JM- and ^(V) = {/GC a (V);ll/IU<«)}. 

s,te[e u e 2 ] r ~ s l 

Using this notation, we define in a natural way 

cfoo = {/Gd(n ii5/iu<oo}, 

and recall that we have also defined a norm || • || Kj oo.t at equation ([7j). In the same way, 
for h G Cs(y), we set 



IWI7.P = SU P 

s,u,te[hM 





hsut 




\u - s\7\t 


— u 


p 



(12) 



WHu = inf jXlll^ll^-^ h = Ylk, < pt < /u| , 

where the last infimum is taken over all sequences {hi, i G N} C ^(V) such that 
h = hi and over all choices of the numbers pi G (0, p). Then is easily seen to be 
a norm on C^{V), and we define 

C$(V) := {heC 3 (V); \\h\\,<oo}. 

Eventually, let Cl + (V) = U^iC^V), and note that the same kind of norms can be 
considered on the spaces ZC^{V), leading to the definition of the spaces ZC^iV) and 
ZC\ + {V). In order to avoid ambiguities, we denote in the following by Af[- ; C K ] the 
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^-Holder norm on the space Cj, for j = 1, 2, 3. For £ £ Cj(^), we also set jV[C; C°(V")] = 
Bup ae[ftA]J ||C*||v. 

The operator 5 can be inverted under some Holder regularity conditions, which is 
essential for the construction of our generalized integrals. 

Theorem 2.1 (The sewing map). Let fi > 1. For any h £ ZC^{V), there exists a 
unique Ah £ C^iV) such that 5(Ah) = h. Furthermore, 

||A%< ^-M[h-C^{V)\. (13) 

This gives rise to a continuous linear map A : ZC^iV) — > C% (V) such that SA = id Z cg(v)- 

Proof. The original proof of this result can be found in [18] . We refer to [TTj IT9] for two 
simplified versions. 

□ 

The sewing map creates a first link between the structures we just introduced and the 
problem of integration of irregular functions: 

Corollary 2.2 (Integration of small increments). For any 1-increment g £ C 2 {V) such 
that Sg £ C\ + , set h = (id — A5)g. Then, there exists f £ Ci(V) such that h = 5f and 

( 5 f)st = tTT lim n X^+i' 

n s t | — , 

where the limit is over any partition U st = {to = s, . . . ,t n = t} of [s, t] whose mesh tends 
to zero. The 1-increment Sf is the indefinite integral of the 1-increment g. 

We also need some product rules for the operator 5. For this recall the following 
convention: for g £ C n {[£ h £ 2 ]R M ) and h £ C m {[£ u £ 2 ]; R d ' p ) let gh be the element of 
C n+TO _i([^ 2 ];R^) defined by 

= 9t lt ...,t (14) 
for ti, . . . , t m+n _i £ [£i, £ 2 ]- With this notation, the following elementary rule holds true: 

Proposition 2.3. Let g £ C 2 ([4, £ 2 ]; WL l ' d ) and h e Ci{[liJ 2 \^ d ) ■ Then gh is an element 
o/C 2 ([£i,£ 2 ];M') and 5(gh) = 5gh-g5h. 

2.2. Random differential equations. One of the main appeals of the algebraic inte- 
gration theory is that differential equations driven by a 7-Holder signal x can be defined 
and solved rather quickly in this setting. In the case of an Holder exponent 7 > 1/3, the 
required structures are just the notion of controlled processes and the Levy area based 
on x. 

Indeed, let us consider an equation of the form 

m 

dy t = a(y t )dx t = J2 a{l) (yt) dx li te[0,T], y = a, (15) 
i=i 

where a is a given initial condition in M. d , x is an element of C7([0,T]; R' m ), and a is a 
smooth enough function from M. d to M. d,m . Then it is natural (see [33] for further expla- 
nations) that the increments of a candidate for a solution to f TTS]) should be controlled 
by the increments of x in the following way: 
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Definition 2.4. Let z be a path in C*(R d ) with 1/3 < k < 7. We say that z is a weakly 
controlled path based on x if Zq = a, with a G R d , and 5z G C 2 (R d ) aas a decomposition 
5z = (5x + r, that is, for any s,t G [0, T], 

(5z) st = ( s (Sx) st + r st , (16) 

with C G C^M ''™) and r G Cf K (R d ). 

The space of weakly controlled paths will be denoted by Q;!: a (R d ), and a process z G 
Q% a {R d ) can be considered in fact as a couple (2, £). The space Q^ a (R d ) is endowed 
with a natural semi-norm given by 

N[z- Ql a (R d )] (17) 

= Af[z; C?(R d )} + JV[C; C°(R d ' m )} + JV[C; C^R^™)] + A/>; C 2 2K (R d )], 

where the quantities J\f[g; Cf\ have been defined in Section 12.11 For the Levy area 
associated to x we assume the following structure: 

Hypothesis 1. The path x : [0, T] — > R m is ^-Holder continuous with ~ < 7 < 1 
and admits a so-called Levy area, that is, a process x 2 G C 2 27 (R m ' m ), which satisfies 
Sx 2 = 5x (8> 5a:, namely 

[(5x 2 ) sil t] (z',j) = [fe*]^^]^, 

/or any s,u,t £ [0, T] and 2, j G {1, . . . , m}. 

To illustrate the idea behind the construction of the generalized integral assume that 
the paths x and z are smooth and also for simplicity that d = m = 1. Then the 
Riemann-Stieltjes integral of z with respect to x is well defined and we have 

z u dx u = z s {x,t — x s ) + / [z u — z s )dx u = z s (5x) s t + / (5z) su dx u 



for £1 < s < t < £ 2 . If z admits the decomposition ffTB"]) we obtain 

t r-t r-t 

(5z) su dx u = / (Cs(5x) su + p su ) dx u = ( s / (5x) su dx u + I p su dx u . (18) 



Moreover, if we set 



:= / (Sx) su dx u , £i<s<t<£ 2 , 

J s 

then it is quickly verified that x 2 is the associated Levy area to x. Hence we can write 

J z u dx u = z s (Sx) sz + (s (x 2 ) st + J Psu dx u . 
Now rewrite this equation as 



/ p su dx u = / z u dx u - z s (5x) st - C (x 2 ) st 

J s J s 



(19) 



and apply the increment operator 5 to both sides of this equation. For smooth paths z 
and x we have 

z dx I =0, 5(z5x) = —5zSx, 
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by Proposition 12.31 Hence, applying these relations to the right hand side of ffl9|) . using 
the decomposition f fl6|) . the properties of the Levy area and again Proposition I2.3[ we 
obtain 

5 I / pdx) = (5z) su (5x) ut + (5() su (x 2 ) ut - C s (<5x 2 ) sut 

\J J J sut 

= Cs(Sx) su (Sx) ut + Psu (Sx) ut + (5()su(X- 2 )ut ~ (s(Sx)su {5x) ut 
= Psu{5x) ut + (5Qsu (x 2 )„t. 

So in summary, we have derived the representation 



5 ( pdx 



Psu 

{5x) ut + (S() 



- sut 



As we are dealing with smooth paths we have 5 (J pdx) G ZC\ + and thus belongs to 
the domain of A due to Proposition 12.11 (Recall that 55 = 0.) Hence, it follows 

/ Psudx u = A st (p5x + 5(x. 2 ) , 

J s 

and inserting this identity into ( Tl8l) we end up with 

J z u dx u = z s {5x) st + (s (x 2 ) st + A st {pSx + 5(x 2 ) . 

Since in addition 

p5x + 5C* 2 = -5(z5x + C* 2 ), 

we can also write this as 

yW(id-M) ( ^ + cx*). 

Thus we have expressed the Riemann-Stieltjes integral of z with respect to x in terms 
of the sewing map A, of the Levy area x 2 and of increments of z resp. x. This can now 
be generalized to the non-smooth case. Note that Corollary 12.21 justifies the use of the 
notion integral. 

In the following, we denote by A* the transposition of a vector resp. matrix, and by 
A\ ■ A 2 = Tr^AiA^) the inner product of two vectors or two matrices A\ and A 2 . 

Proposition 2.5. For fixed | < k < 7 , let x be a path satisfying Hypothesis^ Further- 
more, let z e Q% a ([£i, ^2]; such that the increments of z are given by (TJ7|). Define 
z by Z£ 1 = a with tiel and 

(5z) st = [(id-A5)(z*5x + C-^)] st (20) 

for £1 < s < t < £ 2 - Then J(z* dx) := z is a well-defined element of Q^. & ([£i,£ 2 ];M.) 
and coincides with the usual Riemann integral, whenever z and x are smooth functions. 

Moreover, the Holder norm of J{z* dx) can be estimated in terms of the Holder norm 
of the integrator z. (For this and also for a proof of the above Proposition, see e.g. [18].) 
This allows to use a fixed point argument to obtain the existence of a unique solution 
for rough differential equations. 

Theorem 2.6. For fixed | < K < 7, let x be a path satisfying Hypothesis \J\ and let 
a G C 3 (M a! ; M d ' m ) be bounded with bounded derivatives. Then we have: 
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(1) Equation 173]) admits a unique solution y in Q X K a ([0, T]; M. d ) for any T > 0, and 
there exists a polynomial Pt '■ R 2 — >■ M + sitc/i i/iai 

A% Q:, a ([0,T];M d )] < Pt(|M| 7>00 , t , ||x 2 || 27 ) (21) 

holds. 

(2) Lei F : R d x C7([0,T];R m ) x C 2 27 ([0, T]; R m ' m ) -» C7([0, T]; M d ) 6e the mapping 
defined by 

F (a,x,x 2 ) = y, 

where y is the unique solution of equation (T5\) . This mapping is locally Lips- 
chitz continuous in the following sense: Let x be another driving rough path with 
corresponding Levy area x 2 and a be another initial condition. Moreover denote 
by y the unique solution of the corresponding differential equation. Then, there 
exists an increasing function Kt : K 4 — > M + such that 

(22) 

x (|a — a\ + \\x - x|| 7(00) t + ||x 2 - x 2 || 27 ) 

holds, where we recall that ||/||^ )00j t = ll/lloo + denotes the usual Holder 

norm of a path f G &([0, T\\W). 

Remark 2.7. Inequality l[21\) implies in particular 

\(Sy) st - o-(y s )(5y) st \ <\t- s| 2k P T (||a;|| 7i00iT , ||x 2 || 27 ). (23) 

This estimate will be required in the proof of Lemma \4-3[ 

The above Theorem improves (slightly) the original formulation of the Lipschitz conti- 
nuity of the Ito map F, which can be found in [18J, concerning the control of the solution 
in terms of the driving signal. Therefore (and also for completeness) we provide some 
details of its proof in the appendix. A similar continuity result can be found in [16J, 
where the classical approach of Lyons and Qian to rough differential equations is used. 

2.3. Application to fBm. The application of the rough path theory to an equation 
with a particular driving signal relies on the existence of the Levy area fulfilling Hy- 
pothesis [TJ In our setting, the driving process is given by an m-dimensional fractional 
Brownian motion (B^\ . . . , £?( m )) with Hurst parameter 7 > 1/3. 

To the best of our knowledge, there are three known possibilities to show the existence 
of the associated Levy area B 2 = (B 2 (i, j))jj=i,... ;m : (i) By a piecewise dyadic linear 
interpolation of the paths of B, as done in [B|. (ii) Using Malliavin calculus tools in 
order to define B 2 as a Russo-Vallois iterated integral, similarly to what is done in [30] 
to construct a delayed fractional Levy area, (iii) By means of the analytic approximation 
of B introduced by Unterberger in [37]. Actually, all three methods lead to the same 
Levy area. The equivalence between the first two constructions has been established 
by Coutin and Qian through a representation formula (see Theorem 4 in [6]). The 
convergence results we are going to establish show that the Levy area recently obtained 
by Unterberger in [27J coincide with the previous ones. Note that this question had been 
left open by the author in the latter reference, so that the following Proposition 13.71 has 
an interest in itself (see also [21] for a partial result in this direction). 
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We resort here to the analytic definition of the fractional Levy area, since we use the 
pointwise estimates of [31], which were derived in this setting. Let us recall the main 
features of the analytic approach. 

2.3.1. Definition of the analytic fBm. The article [37J introduces the fractional Brownian 
motion as the real part of the trace on R of an analytic process T (called: analytic 
fractional Brownian motion [35]) defined on the complex upper-half plane Il + = {z £ 
C; ^s(z) > 0}. This is achieved by an explicit series construction: for k > and z £ 

set 

H ^ [MZM l r(2-2H + k) (z + i\ 2H - 2 (z-i\ k 
fk{z) ~ 2 V 2 cos 7rH V T(2-2H)k\ ) \zTi) ' (24) 

where F stands for the usual Gamma function. These functions are well-defined on 
and it can be checked that 

£/*(* + if) 7k (y + if) = K'~ (\ ( m + V2 ) ; x, y 

where K ~ is a positive kernel defined on R+ x 1 x M given by 
K'-(ri;x,y) 

We also set 

K'' + (r);x,y) 

Now define the Gaussian process T' with "time parameter" z £ n + by 

=£/*(*)& (25) 

fc>0 

where (£k)k>o are independent standard complex Gaussian variables, i.e. E[^-^] = 0, 
E[£j£fc] = 5j t k. The Cayley transform z H- ^ maps n + to V, where V stands for the unit 
disk of the complex plane. This allows to prove that the series defining V is a random 
entire series which is analytic on the unit disk and hence the process V is analytic on 
Furthermore, restricting to the horizontal line M. + i|, the following identity holds: 

E[V{x + irj/2)T'(y + i V /2)} = K''~( m x, y). 

One may now integrate the process V over any path 7 : (0, 1) — > U + with endpoints 
7(0) = and 7(1) = z £ n + U R (the result does not depend on the particular path but 
only on the endpoint z). The resulting process, which is denoted by T, is still analytic on 
Furthermore, the real part of the boundary value of V on R is a fractional Brownian 
motion. Another way to look at this is to define T(rj) := {T(t + it]); t £ R} as a regular 
process living on R, and to observe that the real part of r(r/) converges for rj — y to a 
fractional Brownian motion. The following Proposition summarises what has been said 
so far: 

Proposition 2.8 (see [371 135])- Let V be the process defined on U + by relation ( dSp . 



m-2H) 

2cos7ri/ 



■11 



y) + 77) 



2H-2 



H(l-2H) 
2 cos tiH 



(+i(x - y) + rj) 



2H-2 
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(1) Let 7 : (0, 1) — > IT + be a continuous path with endpoints 7(0) = and 7(1) = z, 
and set T z = J T' u du. Then V is an analytic process on Furthermore, as 
z runs along any path in IT + going to t £ R ; the random variables T z converge 
almost surely to a random variable called again T t . 

(2) The family {T t ; t £ IR} defines a centered Gaussian complex-valued process whose 
paths are almost surely K-H6lder continuous for any n < H . Its real part B : = 
{2^RT t ; t £ M} has the same law as fBm. 

(3) The family of centered Gaussian real-valued processes B{rj) := {25Jr t+ir/ ; t £ IR} 
converges a.s. to B in a-Holder norm for any a < H, on any interval [0, T] 
with T > 0. Its infinitesimal covariance kernel 'E[B' x (r])B' y (ri)] is K'(r];x,y) : = 
K'' + (r];x,y) + K' '<- (77; x , y) . 

2.3.2. Definition of the Levy area. Consider now an m-dimensional analytic fBm V = 
(r^, . . . , I^" 1 )). Since the process B{j]) is smooth, one can define the following integrals 
in the Riemann sense for all < s < t < T, 1 < ji, j 2 < m and rj > 0: 

B^( 3l ,j 2 )= / dB^(rj) / dB$\rj). (26) 

J s J s 

It turns out that B 2 ' 77 converges in the Holder spaces C| K from Section |2~T1 (see [371 135]). 
which allows to define the Levy area in the following way: 

Proposition 2.9. Let T > and define ~B? ;n by equation [26]) . Let also < 7 < H. 
Then B satisfies Hypothesis U\ in the following sense: 

(1) The couple (B(r]),B 2 ^) converges in V(Q;€^{[0,T\;R) x C 2 27 ([0, T] 2 ; R m ' m )) for 
all p > 1 to a couple (5,B 2 ) ; where B is a fractional Brownian motion. 

(2) The increment B 2 satisfies the algebraic relation 5B 2 = 5B ® SB. 

One of the advantages of the analytic approach is that an expression for the covariances 
of the Levy area can be easily derived by dominated convergence. We have 

E[B 2 Siti (t,j)B 2 S2t2 ^j)] (27) 

ftl rt 2 PUl PU 2 

= H 2 (2H-1) 2 / / / \ui-u 2 \ 2H - 2 \v l -v 2 \ 2H - 2 dv 2 dv l du 2 du l 

J Sl J s 2 J Si J s 2 

for < S\ < tx < s 2 < t 2 < T and i,j — 1, . . . ,rn. 

Moreover, B(rj) satisfies similar stationarity and scaling properties as the fBm itself. 

Lemma 2.10. We have 

(1) (stationarity) 

{{SB{ V )) S!U+S , 0<u<T-s} = {B{r]) u , 0<u<T-s}, 

(2) (scaling) 

{B(v)cu, 0<u< T/c} = {c H B (J) , < u < T/c} . 

The above Lemma can be shown by straightforward calculations exploiting that B{ji) 
is a Gaussian process with covariance kernel K 1 and will be useful to derive the scaling 
property of the fractional Levy area. See Lemma 13.11 below. 
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3. Approximation of the Levy area 

Let V n .T be the uniform partition {i£ = — , k = 0, . . . , n} of [0, T], and let B n,T be 
the linear interpolation of B based on the points of V n .T- More precisely, B n ' T is defined 
as follows: for t G [0, T], let k G {0, 1, . . . , n — 1} be such that f% < t < t^ +1 . Then we 
have 

^• T = ^+(^f)(^)^ +1 . (28) 

Let also B 2,n ' T be the Levy area of B n ' T , which is simply defined in the Riemann sense 
by 

B^ T (*,j)= f [ ul dB:^dB:^. 

J s J s 

The first step in the convergence analysis of our Milstein type scheme is to determine the 
rate of convergence of the couple (B n,T , B 2 ' n,T ) towards (B, B 2 ). The current section 
is devoted to this step, which can be seen as an extension of [31] to Holder norms. 
Throughout the remainder of this article we will denote unspecified non-negative and 
finite random variables by 8, indicating by indices on which quantities they depend. 
Similarly, we will denote unspecified constants, whose specific value is not relevant, by 
C or K. 

3.1. Preliminary tools. As a first preliminary step, let us state the following elemen- 
tary lemma about the stationarity and scaling properties of the fBm B and its piecewise 
linear interpolation B n,T resp. about the scaling property of the Levy areas B 2 and 

B 2,n,T 

Lemma 3.1. Consider a point s G V n ,T- Then 

{(SB) StU+s ,(5B n ' T ) StU+s ), 0<u<T-s}±{(B u ,BZ> T ), 0<u<T-s}. (29) 
Furthermore, if c > 0, then 

{(B cu , B£f), < u < T/c} k {c H (B u , B^/ c ), 0<u< T/c} . (30) 
Finally, let s,t G V n ,T with s <t. Then we have 

(BUi,j),B 2 s f' T (i,j)) = (t-sr^B^jlBlt^^j)) (31) 
for all i, j = 1, . . . , m. 

Proof. These assertions are of course consequences of the stationarity and scaling prop- 
erties of fBm, i.e. for any c > the process 

5.0 = c H B { J c (32) 

is again a fBm, and for any h G 1R the process 

£0 = {5B^) K+h (33) 

is a fBm. 

Recall that the points of V n ,T are given by tf = — for all j 6 N, and introduce the 
two mappings F™' T and F+ T defined on R + by F™' T (n) = t™ and F™ ,T (u) = tf +1 if 
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tf < u < tf +1 . With these notations, one has B™' T = G n,T (B) u , where the measurable 
mapping G n ' T : C(R + ; R m ) -> C(R + ; M m ) is defined by 

u - F?' T (u) ( \ 
G n > (y) u = V F n,T {u) + ~ y F-- T (u) ) » m G K + . 

Now, in order to establish f l29l) . note that F± T (u + s) = F± T (u) + s if s G "P n ,T- It is 
then easily seen that 

((6B) s> . +s ,(SB n > T ) Sj . +s ) = ((5i?) s , +s ,G"' T ((5 J B) s ,. +s ), 

so that, due to the stationarity property of fBm, the following identity in law for processes 
holds true: 

((8B) S ,. +S , (8B n > T ) s , +s ) k (£>, G n ' T (B)) = (B,B n ' T ). 

The proof of f)30p is quite similar. In fact, one has F± T (c ■ u) = c ■ F± T ^ c (u) and so 
B™^ = G n ' T ' c (B c .) u . Thus it holds, thanks to the scaling property of fBm, 

(B C .,B:: T ) = {B c .,G n < T / c {B c .)) k (c H B,G n ' T / c (c H B)). 

Identity fl30l) is then a consequence of the linearity of G n,T l c . 

Now it remains to establish (I3~TI) . Note first that Proposition 12.91 implies that 

(BLB^)=lim(B(^,B(,);r T ) (34) 

in probability. Here B{j]f'^ ,T is the Levy area associated to the piecewise linear inter- 
polation of B(rj) with stepsize T/n. 

Since B{rj) is analytic, the above Levy areas can be approximated by a standard Euler 
quadrature rule, i.e. we have 

B( V ) 2 st = lim l k (B( V )l) B( V ) 2 f' T = lim l k (B( V ) 2 f' T ) (35) 

fc— >oo fe— >oo 

almost surely, where 

k 

UB(v)l) = £ {( SB (v))s, i{t - s)+s } ® {(SB(r))) Ut _ s)+smt _ s)+s ) 

i=0 

Using again the G n,T notation and setting r] st = 73- we have 

{UB^DMBitff- T )) = 

k 

(E{ww}®{ 

i=0 

fe 

E GB,T ((^('?))...+.)i(*-.) ® {c n ' T ((^('7))^)^(t- s ) -^ T ((^(^)) Si . +s ) |(t _ s) }). 

i=0 



L -(t-s)+s 



(5B(rj)) 
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Thus, invoking Lemma [2.101 and setting t] st = 7^-, we end up with 

k 
i=0 

fc 
i=0 

fc 

S^(^)i(t- S ) ® (^))|(t- s ),i±I(t- s ), 
i=0 

fc 

i=0 

it 

fc fe ' fe 

(t _ S )^^G"' T /^ (5 ( v st )) ± <g> (^ T /('- s) (5 (77 s *))). 



i=0 



i=0 



that is 



(^(B^X^)^) 



Clearly, we also have 



(t-s) 2H [T k [B 



B 



£ — sJ 01 



,X fc ( B 



£ — s/ 01 



2,n,T/(t-s) 



V 



B 



£ — S/ 01 fc-s-oo 
^ n 2,n,T/ (t-s) 

t — S ) 01 fc-s-oo 



lim Xfc ( B 
lim Xfc ( B 



V 



£ — s/ 01 

77 \ 2,n,T/(t-s) 



almost surely and 



R 2 R 2,n,T/{f-s) 

■ D 01' - D oi 



) = i im . ( B 



£ — s/01 

„ \ 2,n,T/(t-s) 

B' 



r?->0 V \£ — s/01' \£ — s/01 
in probability So, combining ( 134)) . (155]) . ( |36|) . (137|) and ( 138|) . we obtain 



E 



V ( B si , B^ 



lim E 

fe—> 00,77— >o 



^(X^B^^^^B^)^)) 
\p({t - S ) 2 ^X,(B (^)o0> (* " s) 2H lk(B ( V st ) 

v ((t- s rB 2 01 ,(t-sr H B 2 t T/it - s) 



lim E 

fc— > 00,77— >o 



E 



2,n,T/(t 
01 



for any function ip G Cb((M. m (g) M" 1 ) 2 ), which concludes the proof of ( 13~TT) . 
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The next auxiliary result is an upper bound of the modulus of continuity of fBm and 
is a consequence of Theorem 3.1 in |38j . 

Lemma 3.2. Let T > 0. There exists h* > and a finite and non-negative random 
variable 0H,h*,T such that 

sup \(5B) t , t+h \ < 6 Hih *, T -h H -y/\ log(l//i)| 

te[0,T-h] 

for all h £ (0,h*). 

The classical Garsia lemma reads as follows: 
Lemma 3.3. For all 7 > and p > 1 there exists a constant C 7 . p j > such that 

( r T r T usf) \ 2p \ 1/(2p) 

AT[/ ; C7([0,T];lR')]<C 7 , p ^y o jf i^j^ ^d.j 

/or aH/edl^T];!'). 

Finally, we also need to control the Holder smoothness of elements of C2, beyond the 
case of increments of functions in C\. The following is a generalization of the Garsia- 
Rodemich-Rumsey lemma above. 

Lemma 3.4. Let n > and p > 1. Let R £ C 2 ([0, T]; M.^ z ) wit/i £ C* ([0, T]; R 1 -'). // 

dudv < 00, 



'0 ^0 





Ruv 


\2p 


\u 


— v\ 


2np+2 



then R £ C$ ([0, T]; R ' ). In particular, there exists a constant C K)P j > 0, swc/i t/iat 



/ f T f T ID |2p \V( 2 P) 





Ruv 


\ 2P 




— v\ 


2np+2 



3.2. Approximation results. Recall that our aim here is to show the convergence of 
the couple (B n ' T , B 2,n ' T ) towards (S, B 2 ) in some suitable Holder spaces. A similar 
result was obtained in [6], but with the following differences: (i) The authors in [6] 
studied the p-variation norm of B 2 — B 2,2 "' T using dyadic discretisations, while we are 
working in the Holder setting, (ii) The rate of convergence for the approximation was 
not their main concern, and the convergence rate stated in [6j Corollary 20] is not sharp. 

Let us now start with a first moment estimate for the difference B 2 — B 2,n,T , for 
which we will use the error bound for a trapezoidal approximation of B 2 derived in |31j . 
Moreover, recall that we denote by B n,T the piecewise linear interpolation of B on [0, T] 
with respect to the uniform partition V n< T = {£%) k = 0, . . . , n}, where t\ — — , and by 
B 2, ™' T the corresponding Levy area. 

Proposition 3.5. Let p > 1 and H > 1/4. Then, we have 

p\ Vp 



E 



p. 2 p>2,n,T 
■ D 0,T - D 0,T 



<K n -T 2H -n- 2H+l l 2 . 
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Proof. First note that the random variable Bq T — Bq^' t belongs to the sum of the first 
and the second chaos of B (we refer to [32] for a specific description of these notions). 
So all moments of Bq t — Bq'^' t are equivalent and it suffices to show that there exists 
a constant K > such that, for all i,j — 1, . . . , m, 



E 



2.n,T 



B o,r B 0) , 



1/2 



(39) 



Consider first the diagonal elements of Bq t — Bq'^' t . In this case, we have Bq T (j, j) 
{Bff/2 and 

/" T B n,T,(j) rfB?™ 
Jo 

n— 1 ra— 1 ^>t™ 

l k l k l k+l i-^d 



B2,n,T/ ■ 
or CM; 



fc=0 
n-l 



k=0 Jt k v 7 



S ^2*2+! + 2 ) = 2 

fc=0 v 7 



Hence it follows 



I B®dB®-f Bf(^BfW=0. (40) 
Jo Jo 



Now consider the off-diagonal terms of Bq T — Bq'^ t . Without loss of generality we 
can assume that i > j. Proceeding as above we have 



" " 2 ^ ' fc L k L k + l 

k=0 



Thus, [3U Theorem 1.2] can be applied and yields 



E 



b o,t(m) -Bj'? ,T (i,j 



1/2 



< K ■ T 2H ■ n~ 2H+1 ' 2 . 



(41) 
□ 



The next result gives an error bound for the piecewise linear interpolation of B. Note 
that similar estimates as in the next lemma can be found in [10], where the case H > 1/2 
is considered. 

Lemma 3.6. Let < 7 < H. Then, there exists a finite and non-negative random 
variable 0h,^,t such that 

M[B n > T - B ; C7([0, T])] < H „? ■ v/bgM • 

for n > 1 . 

Proof. Clearly, we have to find appropriate bounds for 

\5(B n > T - B) st \, s,te[0,T}. 
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First note that there exists a strictly positive xn n such that the mapping / : (0, T] — > 
[0, oo), f(x) = x"-V|log(l/x)| 

is increasing on x G (0, Xff, 7 ). Without loss of general- 
ity, we assume that T/n < inf (a;# i7 , h*), where h* is defined by Lemma [3.21 

(i) First, consider the case where \t — s\ > — . Let us assume also without loss of 
generality that i£ < s < t% +1 < tf < t < tf +1 for some k < I and recall that f% = kT/n. 
Then 

B^ T = B t » + (^*) 5B t?J? _ and B? T = B t? + ^—^ I 5B t?t?+i , 



T/n J k k+1 1 1 V T / n 



so that 

>n,T 



ISiB 11 ' 1 - B) at \ < \5B t n s \ + \5B tl J + \8B t n t n +i \ + \5 B tftf J 

< 40tf iT V|log(7i/T)| (^j < 6 HtT \t - s|V|log(n)|n-( H - 



7) 



using Lemma [3. 2 [ 

(ii) Now, suppose that \t — s\ < T/n with for instance < s < t < t% +1 . In this case, 
and thus 

}n,T m I ^ irn I , |xr?n,T| 



H-l 



15(5™ - B) st \ < \5B st \ + \8B st 

< e HtTy /\log(l/(t-s))\\t - s\ H + 9 HtT \t - s|V|log(n/T)| 

< e HtTy /\\og(i/(t-s))\\t - s\ H + e H>T \t - s |VliogHI^ (H-7) - 

Using the monotonicity of x H- x 11 ^ ^/| log(l/rr)|, it follows 

\5(B n > T - B) st \ < 9 H>T \t - sP^fhgJnj\n-^\ 

(iii) The same estimate as above also holds true if \t — s\ < T/n and f% < s < t'^ +1 < 
t <■ t n 

(iv) Combining (i)-(iii) yields the assertion. 



□ 

Now we determine the error for the approximation of the Levy area. 

Lemma 3.7. Let 1/4 < 7 < H. Then, there exists a finite and non-negative random 
variable 9h,^,t such that 



^ B 2,n,r_ B 2. c 27 ([0;T]) ] < Q H ^ T . y / bg(n) • n~^ H ~^ 

for n > 1 . 

Proof. In this proof we will denote constants (which depend only on p, q, e, 7 and T) 
by K, regardless of their value. 

Step 1. We will first show that 

(E |A/"[B 2 ^ t - B 2 ;C 2 2 \[0,T])]\ q y /Q < K ■ (n~ 2 ^ + n~ H ). (42) 
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For this, we have to consider the family of increments A n ' T (i, j) G C%, defined by 
A n /(iJ)= f{8B®) M dB®- f {5B n ™) su dB n u T ^ 

J s J s 

for i, j = 1, . . . , m. By symmetry we can assume 1 < j < % < m. 

We distinguish several cases for s,t G [0, T]. 

(i) Assume that \t — s\ > — and s,t G V n) T, i-e. s = — and t = — for k < I. Then the 
scaling properties of fBm, see Lemma I3.1[ yield 



t-s 



u u 



t-s 



» dB^ 

jn,T/(t-s),(i) dB n,T/(t-s),(j) \ _ 



Since -r— = -rr we have 

t—s l—k 

| B n,T/(t- s ),(i) j u e [0> x] | = n G [0,1]}. 

Now Proposition 13.51 gives 

(E \A:i T (i,j)\y /P <K-\t- s\ 2H ■ \l - kl 2H+1 ' 2 <K-\t- s^ 2 ■ n~ 2H+1 ' 2 

< K ■ \t — s| 27 • rT 2 ^ H ~ 1 \ (43) 

with 7 > 1/4. 

(ii) Assume now that (t — s) > — with s < t% +1 <tf <t < tf +1 . Using the cohomologic 



relation 5(5A n ' T (i, j))st» +1 t»t = 



^ n si T ihj) = A n J + A^f t Ji,j) + A$(i,j) 

I I 

\n,Ti 



0, we obtain 



+ 5{A^ (i,j)) sW + 5(A n '+ (i,j))tz +1 t T t. (44) 



For the term A^^ t n(i, j), we can use the first step to deduce 



E 



p\ i/p 



A'^rU.j) ) <K.\tf- tlX ■ ^ 2(H ^ <K-\t-sf<. n~ 2 ^\ 



1+n 



To deal with the last two terms of (144]) . remember the algebraic relation 

5(A n > T (i, j)) = 5B {i) ■ 5B U) - 5B n < T > {i) ■ 5B n > T > {j \ (45) 

which entails here 

m n ' T ^j)) sq+lt \ < |(*BW)rfjJ ■ \(6B«% +it \ + \(5B^) sq+i \ ■ \(5B^) n+it l 
and we easily get 

(E \5(A n ' T (iJ)) stt+it \ P ) 1/P <K-\t-s\ H - (T/n) H < K • \t — s| 27 • (n~ 2 ^ + n~ H ). 
Similarly we obtain the same estimate for E[|<5(A ra,T (z, j))t™ +1 t] i t\ p ] 1 ^ p - 
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As for the term A n s [ n one has, on the one hand 



E 



■fc+i 



(6B®) su dBW 



\t n — <? 

I L k+1 b 



2H 



E 



u u 



P\ 1/p 



<K-\t-s\^-n- 2 ^\ 



(46) 
(47) 



where 7 < H. On the other hand, 



'fc+1 



SB*?'® dB*™ 



SB n,T,(i) SB n,T,(j) 

L k L k °k+l 



(T/n)> 

So for 7 < H, an application of the Cauchy-Schwarz inequality yields 



5B n,T,{i) SB n,T,U) 

U J-J j-n. -Mil V ±Jj.rij.n 

L k b k+l L k L k + 1 



E 



•fc+i 



{5B n ™) su dB% 



T,(j) 



P\ l/i> 



< A" • |i - s\ 



27 . n -2(H- 7 ) 



(48) 



Putting together relation flU} and (gBJ), we obtain (E[|^J ^z, j)I p ]) Vp < ^ |t - s| 27 
n -2(H--y)^ Furthermore, the term A^(i, j) can be handled along the same lines, 
(iii) It only remains to analyze the case (t — s) < — . For i£ < s < t < t^ +1 we have 



E 



and 



and thus 



(6BM) m dB® 



5B n T l{i) dB n T,(j) 



< K • \t — s\ 2H < K • \t — s]' 21 • n 



2 7 -2{H-i) 



< 



(t-s) 
2{T/n) 



X B n,T,(i) 

UJ-Jj.nj.rt 

L k L fc+1 



5B 



n,T,(j) 

j.n .in 
L k L k+l 



E 



{8B n ™) m dB^ 



■ (./) 



P\ 1 /p 



< A- \t-s\ 21 -n 



27 . _-2(H-tO 



The case (t — s) < — and < s < t% +1 < t < t^ +2 can be treated analogously, 
(iv) Combining steps (i)-(iii) yields that 

(E \AF(i,j)\y /9 <K-\t- s| 27 • (n~ 2 ^ + n~ H ) 
for all s, t G [0, T] and 1/4 < 7 < A". 



(49) 



Step 2. Before we can apply Lemma [3 A\ we need additional preparations. First, notice 
that (I45p can also be written as 

5(B 2 - B 2 ' n > T ) =[5(B- B n ' T )] ®5B + 5B n ' T ® [5 (B - B n ' T )] , 

so that 

m 2_ B 2,n,T U] 

< \t - u\~<\s - u\i (2N[5B; C\\ ■ M[5(B - B n < T ); C 7 ] + (j\f[5(B - B n > T ); C 7 ]) 2 ) 

and thus 

Af[5(B 2 - B 2 '"< T ); C 27 ] < 2Af[5B; C\\ ■ M[8{B - B n ' T ); C 7 ] + (M[6(B - B n > T ); C 7 ]) 2 . 
Lemma [3.61 now gives 

-(H-7) 



AA[5(B 2 -B 2 ^);C 27 ] < 0H l7tT ■ Vtogi 



n) ■ n 



(50) 
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Step 3. Using ([3D]) . Lemma [3T41 entails 
A/"[(B 2 -B 2 '"' T );C 2 27 ([0,T]] 

/ pT pT I /p.2 _ p>2,ra,T\ \2p \ VO) 

- K \L I \u-v\J? dudV ) +K.6 w .^n).n^\ 
for all p > 1. To finish the proof, it remains to show that 



\R n , P \ < 9 7>H>T ■ Vlog(n) • (51) 



where 

: T r T |(b 2 -b 2 ^u 2 ? n 1/(2p) 

it ni p = I / / ; du dv 



o Jo 



it - v\^p+ 2 



However, using ( )49l) with 7 + e/2 instead of 7, we have 
FIR \*< f T /- T E|(B 2 -B 2 ^U| 2 ^ 



</o 



< K / |M \ ^ 9 du cfe • ( n -4(H- 7 - £ /2) P + n -2ffp) 
" 7o Jo \u-v\^p+i { h 



i.e. 



(E|i? n , P | 2p ) 1/(2p) < X / / \u-v\ 2p£ - 2 dudv(n- 2 ^ H -^ + n- H ). 
Jo Jo 

So for p > -, it holds 

(Eli^l*) 1 /^) < # • (n~ 2 ^^ +£ + n~ H ). 

Now, set a = min{2(iJ — 7) — e,H} and let 5 > 0. From the Chebyshev-Markov 
inequality it follows 

F,l R l 2 P „-2pe 

P(n-K,,| > 5) < ^g^^(-) < 
Since p > 1/e we have 

00 

^PK- £ |i? niP |>(5)<oo 

n=l 

for all <5 > 0. The Borel-Cantelli Lemma implies now that n Q_e |-R niP | — > a.s. for 
n — > 00, which gives floTj) by choosing e > appropriately, since 

a — e = min{2(i7 — 7 — e), H — e} > H — 7. 

□ 

Recall that the Wong-Zakai approximation Z n of Y has been defined at equation (jSJ) 

by 

Z'; = a + Y, v {t) (K) dB^ T , t G [0, T] , a G M d . (52) 
i=i ^ 
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In particular, Z can be expressed as Z = F(a, B n,T , B 2 > n > T ) , using Theorem 12. 6 1 Hence, 
as a direct application of Lemmata 13.61 and 13.71 and invoking the Lipschitzness of F, we 
obtain the following error bound for the Wong-Zakai approximation. 

Proposition 3.8. Let T > and 1/3 < 7 < H. Then, there exists a finite random 
variable rj^ a T such that 

\\Y - Z n \\^ T < ■ v/bgM • n-<*-rt 

for n > 1 . 



4. DlSCRETISING THE WONG-ZAKAI APPROXIMATION 

In the last section we have established an error bound for the Wong-Zakai approxima- 
tion Z of the real solution Y . As mentioned in the introduction, the Milstein scheme 
corresponding to Z is exactly our simplified Milstein scheme (JSJ). Thus, it remains to 
determine the discretisation error for Z itself. To this aim, we first give a general error 
bound for the Milstein scheme for ordinary differential equations (ODEs) driven by a 
smooth path x. Since Theorem 12.61 allows to derive a non-classical stability result (in 
7-Holder norm) for the flow of an ODE driven by a smooth path, we can follow here the 
techniques of the numerical analysis for classical ODEs. In a second step, we will apply 
these bounds to our particular fBm approximation. 

4.1. The Milstein scheme for ODEs driven by smooth paths. In this section, 
consider a piecewise differentiable path x £ C([0, T] ; M. ) and a function g e C 3 (M. d ; M. d ' 1 ) 
which is bounded with bounded derivatives. For the ordinary differential equation 

1 

y t = J29 { ' l) (yt)dxf ) , te[0,T], aeR d , (53) 

i=l 

the classical second order Taylor scheme with stepsize T/n reads and 




where = Ylp=i 9p^ p , and where we have set z% = z^n with t% = kT/n. For notational 
simplicity we will write in the following tk instead of f%. Introducing the numerical flow 




we can write this scheme as 

z£ = a, zl +i =ty(z%;t k ,t k+1 ), fc = 0,...,ra-l. 
For q > k we also define 



t k , t q ) := (■; t,_ l9 t q ) o (•; t q _ 2 , ■ " • ^{z; t k , t k+1 ). 
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Moreover, the flow s,t) of the ODE (153]) is given by s,t) := y t , where y is 
the unique solution of 

i 

y t = Y,9 {i) (yt)dxf\ te[s,T\, y s = z. (56) 

A straightforward Taylor expansion of the flow of the ODE gives that the one-step 
error 

r k = t k , t k +i) - ^(z; t k , t k+1 ) 

satisfies 



<C- sup ||P«I?^V P) l|oo-M^ +i 

i ,j,p=0,...,m 



(57) 



with 



M st 



Furthermore, considering the smooth path x as a rough path, Theorem 12.61 directly 
yields the following stability result for the flow: 

Proposition 4.1. Let 1/3 < 7 < 1 and set ||x|| 7 = ||x|| 7 + ||x 2 ||2 7 . Then, there exists 
an increasing function Ct '■ K — > R+ such that 

\($(z;s,t)-$(z;s,t))-(z-z) 



t-sP 



< C T (\\x\\,) ■ \z - z\ 



and 



\$(z; s,t) - $(z; s,t)\ < Cr(||x|| 7 ) ■ \z — z\ 
for all s, t £ [0, T] and z,z £ R d . 



(58) 
(59) 



The following stability result is crucial to derive the announced error bound for the 
Milstein scheme. 

Proposition 4.2. Let x £ C([0,T];M. 1 ) be a piecewise differentiable path, and g £ 
Cfr(M. d ; M. d ' 1 ). Consider the flow $ given by equation l[5b}) and the numerical flow \1/ 
defined by relation 153]) . For k — 0, . . . , n, let t k — kT/n, y tk = $(a; 0,t k ) and z tk = 
\I/(a; 0,t k ). Moreover recall that we have set 

3 



M\ 



St 



< s < t < T. 



Then, there exists an increasing function Ct '■ R — > K + such that we have 

q-l 



\y tq -z-\ < <?t(||x|| 7 ) • 5>* tfe+i 



\S(y-z n ) tptq 

for < p < q < n. 



k=0 

< <5r(||x|| 7 ) ■ + \t q ~ t P V ■ 



(60) 
(61) 
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Proof. We will use the classical decomposition of the error in terms of the exact and the 
numerical flow: Since z k = §(z k ;t k ,t k ) and y tk = $(zq ; t , one has 

9-1 

y tq - z n q = <S>(z%-t Q ,t q ) - <f>{z%;t q ,t q ) = £ ($(z£;t fc ,t,) - $(4 +1 ;t fe+1 ,t 9 )). 

k=0 

Furthermore, thanks to the relation 



$(z£;t fc ,t g ) = $($(z^;t fc ,tfc + i);t fc+ i,t,), 
the stability result fl59|) implies 

- $(*fc+i;*fc+i,*?)| < Ct(||x|| 7 ) • \$(z%;t k ,t k+1 ) -zl +l \ . 

However, (1ST)) gives 

|d>(^; t fc , - z n k+1 1 = |$(*»; t k , t k+x ) - (z»; t fc , | < C7 • M* tfc+i , 

from which (16T))) is easily deduced. 
Moreover, for q > p we also have 

% - z n )t„t, = t P , t 9 ) - yt P ) - (*(^ n ; t p , t q ) - z n p ) 

= {Hvt P ;t p ,t q )- ytp ) - ($(z;;t p ,t q ) -z n p )) - (y( z ;-,t p ,t q )-<s>(z;-,t p ,t q )) . 

Analogously to the derivation of (1ST))) , one can show that 

M$M) ~ < C ■ Or(||x|| 7 ) ■ J2 M t k t k+1 - (62) 

fc=p 

Using (1551) and fl60|) we trivially end up with (1ST)) . 

□ 



4.2. Application to fBm. In order to apply Proposition 14.21 to the Wong-Zakai ap- 
proximation Z given by ( 152)) note once again that our Milstein-type scheme Z™ = a 
and 

Z? = Z n t + V a w [Z n t ) 5B& + - V V>(Z t n ) 5R ( 2 
«=i m=i 

is obtained by discretising the Wong-Zakai approximation with the standard second 
order Taylor scheme with stepsize T/n given by (1541) . In fact, doing so we obtain the 
numerical flow 



9(z; t k , t k+1 ) :=z + J2 ^{z^C + £ V ^ U) ^ / dB ^ 



Since B n ' is the piecewise linear interpolation of B on [0,T] with stepsize T/n, the 
above iterated integrals can be now expressed as products of increments of B. Indeed, 
according to the fact that 

*Bg? T = 5B^?0, K< T = ^(6B) tktk+1 for ue(t k ,t k+1 ) : (63) 
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it is readily checked that 



8Bf)£ T = 5B® , and / ' ^ SB^ n ' T dB^ T = - SB® 5B® . 



' I I: 

Moreover, invoking relation (163]) and Lemma [3 .2[ we get 

<"tk+l 



\B^ T \ du 



<B H , T n- H [login)] 1 / 2 



for n large enough. Consequently, relation f|6Tj) yields 

sup '^""f!^ 1 < W ^ 3 " +1 [log(n)] 3 / 2 (64) 

p,q=0,l,...,n—l,pj^q |Cp £g | 

for all 7 < if and all n large enough. 
This gives in particular 

sup l6{ ^ Zfl tptql ^ 6 "^t [log(n)] ^ (65) 

p,q=0,l,...,n-l,p^q \tp tq\ 

for 1/3 < 7 < H. 

Now it remains to "lift" this error estimate to [0, T]. For this we need the following 
smoothness result for the Wong-Zakai approximation. 

Lemma 4.3. Let T > and recall that Z™ is defined by equation ( T5l|) . Then there 
exists h** > and a finite and non-negative random variable 0n,h**,a,T such that for all 
h G (0, h**) and all n > we have 



sup \{ST) t>t+h \ < 6 Hjh ~ atT ■ h H • V|log(l/fc)|. 
te[0,T-h] 



Proof. As already mentioned in the proof of Lemma 13.61 note that there exists xh > 
such that the map x n- x H \/| log(l/x) | is increasing on (0,x#]. Set h** = min(x#, /i*), 
where h* is defined by Lemma [3. 2[ and let s,t £ [0,T] such that \t — s\ < h**. 
(i) From ( 1B^|) and (I2"3"j) . we deduce 



K^^-a^)^^'),,! < |t- S r K G(||B^ 

for l/3<K<7<if and an increasing function G : 1R — > M + . Choosing n, 7 sufficiently 
large, we obtain 



\{ST) st -a{X){SB n > T ) st \ < e HM , a>T \t-s\ H ^jlog(^j^ 

(ii) Assume that t L < s < t < t t+1 . One has 

\o-(Z n s )(5B n > T ) st \ < e Hjh *, a , T -\t-s\- {n/Tf- H ^\\og{n/T)\. 
Since \t — s\ < T/n, i.e. n/T < l/(t — s), it follows 

\(ST) st \ < Qh,h* ,&,T ■ (t - s) H • v/|log(l/(t-s))|. 

(iii) Now let < s < ti < t p < t < t p+ i with I < p. Then 

(SB"> T ) st = (B?> T - B tp ) + (SB) tltp + (B tl - B:> T ). (66) 
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As in the proof of Lemma I3.6[ this easily yields 

\a(Z n s )(5B n > T ) st \ < 0* ihv , T • (t - s) H • v/|log(l/(t- S ))| (67) 
for \t — s\ <T/n. Whenever \t — s\ > T/n, decomposition ( 166]) gives 



\a(Z n 8 )(5B n < T ) 8t \ <2^, CT , T -(T/n) ff v /|log(n/T)|+^ iCTiT -(t p -tO H Y|log(l/(t p -t i ))l 

Using that x (-)■ x E \/\ log(l/x)| is increasing, relation ( 1671) is easily recovered, 
(iv) Combining the steps (i)-(iii) yields the assertion. 



□ 



Proposition 4.4. Let T > and 1/3 < 7 < H . Then, there exists a finite and non- 

ble rf^ 1 a j< such that 

Z n - Z n \\^ T < 4\ aT ■ v/Io^H • 



(2) 

negative random variable rj Hn aT such that 



for n > 1 . 

Proof. Denote by U n the piecewise linear interpolation with stepsize T/n of the Wong- 
Zakai approximation Z . Proceeding as in the proof of Lemma 13.61 and using Lemma 
14.31 we have 

\\u n - z n \\^ T < e H „,„,T ■ v/iogR • n-^-A 

Thus, it remains to consider the difference between U n and Z n . For t G [t^, tfc+i] for 
some k we have 

U? - Z? = X - Z? + t - 7 ^5 (Z n - Z n ) t t . 
Assuming additionally that s G [ti,ti + i] and t G [tfc, tfe+i] for some I < k, we have 
5(U» - Z n ) st = 5(T - Z n ) tltk + - Z n ) tktk+1 - - Z n ) tltl+1 . (68) 

(i) Assume that I + 1 < k. Applying f|65l) to relation ( 168]) and according to the fact that 
(s — ti) < T/n, (t — tk) < T/n, we obtain 

\5(U n - Z n ) st \ < 9 H „, aiT \t - sp ■ rT^^y/^nj. (69) 

(ii) Assume that I = k. Here (|68|) simplifies to 

5(£T - Z n )st = - Zn )^ k+1 

and thus ( 1651) combined with the fact that |t— s\ <T/n gives an estimate of the form ( 169] ) 
again. 

Finally, the case = I + 1 can be treated in a similar manner, and this completes the 
proof. 

□ 

Remark 4.5. Putting together Propositions 13.81 and 14.41 our Main Theorem 11.11 now 
follows. 
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4.3. Optimality of the error bound. Reviewing the steps of the derivation of our 
main result, one realises that the final convergence rate n~^ H ~^ ylog(n) is directly 
linked to the error (measured in the 7- Holder norm) of the piecewise linear interpolation 
of fractional Brownian motion. All other estimates lead to higher rates of convergence. 
As a result, in order to prove the optimality of our result, it is natural to consider the 
most simple equation 

dY t {1) = dB[ X \ t G [0, T] , Y = aeR, 

for which our Milstein-type approximation is given by Z n = B n,T . 
First, observe that 



\Y - Z n \L^ T = \\B^ - B (1) ' re ' T L,oc,T > sup 



s,te[o,T] 



\S(B^ - B^ T ) st \ 

\t - S|T 



> 



sup 



| B (D _ B (D,n,T| 



te[o,T] t 1 

Using the scaling and stationarity properties of fBm, we get 

ID n n ' T \ ID D n il I ID D"," 

sup = sup 1 — = 1 ' sup n 



t£[o,T] t 1 te[o,i] T~*t~< te[o,n] n 7 ^ 7 

> n^ H -^T H ~^ sup \B t - £?H = „-(«-7) T h-7 sup \B t -Br 1,n '\ (70) 

te[l,n] te[0,»-l] 

Now let us recall the following result of 



Vn I Id d'M I 1 £ > / 

sup \B t -B t ' - o n v n I — y G, 



<Tr, 



*e[o,i] 



where G is a Gumbel distribution, lim^oo , = = 1 and lim n ^ 00 n H cr n = c#. This 

V 21 °s( n ) 

implies in particular 

TJ 

sup \B t — -B™' 1 ! ——y V2ch- 



y/^og(n) te[o,i] 
Applying again the scaling property of fBm gives 

1= sup \Bt - B? n \ A Va 

log(n) ie[o,n] 

and so 



ch 



\/5og(n) te[o,n-i] 
Going back to (17D1 . this finally yields 



sup 



lim P{£(n) ■ \\Y - Z n || 7i00jT < 00) = 0, 

n— >oo 

if 

liminf £(n) ■ \/log(n) ■ n~^ H ~^ = 00, 
which corresponds to our claim at Remark 11.61 
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5. Numerical Examples 

In the introduction, we stated the conjecture that the error in the supremum norm of 
our proposed modified Milstein scheme satisfies 

\\Y - Z n \\^ T < VH , a , T ■ y/bgfi • (n' H + n~ 2H+1 l 2 ) . 



* mod milstein scheme 
conv. order 2H-1/2 



* mod milstein scheme 
conv. order 2H-1/2 



10 10 10 10 10 10 



10 10 10 10 10 10 



* mod milstein scheme 
conv. order 2H-1. ; 2 



* mod milstein scheme 
conv. order 2H-1/2 



Irs - 



10 10 10 10 10 



10 10 10 10 10 10 



Figure 1. Equation (1711) : pathwise maximum error vs. step size for four 
sample paths for H = 0.4. 



Note that if U n denotes the piecewise linear interpolation of Z with stepsize T/n, then 
we have 

\\Y - U n \\^ T < Tjjj a r p ■ Vlog(n) • n~ H , 

which follows from a straightforward modification of the Lemmata 13.61 and 14.31 Since 
furthermore 

\\Y-Z n \\ 00 , T <\\Y-U n \\ 00 , T + max \Y kT/n - Z% T/n \, 

k=0,...,n ' 

it suffices to consider the maximal error in the discretisation points, i.e. 



max \Y kT/n - Z% T/n \, 

k=0,...,n 



to support our conjecture. 
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^1 ' ' , I Ifl-'l I 

10 to" 1 10 J 10 J 10 ' 10* 10 10 J 10 J 10 ' 10 ' 10' 

stepsize stepsize 



FIGURE 2. Equation ( I7TT) : pathwise maximum error vs. step size for four 
sample paths for H = 0.7. 

Our first example will be the SDE 

dY t = cos(Y t ) dB { t 1] + sin(Yt) dBf\ t £ [0, 1], Y = l. (71) 
Figure 1 shows the maximum error in the discretization points, i.e. 

max^ \Y kT/n (uj) - Z% T/n (uj)\, 

which for brevity we call in the following maximum error, versus the step size 1/n for 
four different sample paths u £ Q for H = 0.4, while Figure 2 shows the maximum 
error versus the step size 1/n for four different sample paths u £ Q for H = 0.7. (So 
small values on the x-axis correspond to small stepsizes, while small values on the y-axis 
correspond to small errors and vice versa.) 

The numerical reference solution is obtained by using our Milstein-type scheme with 
very small stepsize. Since we use log-log-coordinates, the straight lines correspond to 
the convergence order 2H — 1/2. The stars correspond to the error of the Milstein-type 
scheme. For H = 0.4 the estimated convergence rates are in acceptable accordance with 
our conjecture, while for H = 0.7 they are in good accordance. 

As second example we consider the linear equation 

dY^ = Y} 2) dBf\ dY t (2) = Y t (1) dB?\ t £ [0,1], F (1) = 1, F (2) = 2. (72) 
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FIGURE 3. Equation (1721) : pathwise maximum error vs. step size for four 
sample paths for H = 0.4. 

Figures 3 and 4 show again the maximum error versus the step size for four different sam- 
ple paths for H = 0.4 and H = 0.7, respectively. Again the estimated covergence rates 
are in acceptable accordance with our conjecture for H = 0.4 and in good accordance 
for H = 0.7. 

Note that the convergence order 2H — 1/2 is quite slow for small H. In particular, 
for H = 0.4 the convergence order equals 0.3. We suppose that this effect also causes 
the fluctuating behaviour in the estimated convergence rates in the case H = 0.4. 

6. Appendix: Proof of Theorem 12.61 

6.1. Existence and uniqueness of the solution. This section gives some details of 
the proof of point (1) of Theorem 12.61 in the case 7 < 1/2. The case 7 > 1/2 is simpler 
and thus omitted. 

The solution to equation (j!5p is obtained via a fixed-point argument, which is first 
applied locally and then extended to the whole interval [0,T]. 

Notations. For Q X R a ([£i, £2}', K d ) we will write in the following only Q*{[£i, £2]) to simplify 
the notation. In particular, note that the norm Af[-; Q« ([-^i,-^2])] does not depend on 
a G M, d . Moreover, for y G Q%([£i,£2\), which admits the decomposition 

(Sy) st = C a ($£)st + r s t, 
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mod milslein schems 
-conv. order 2H-1/2 



mod milslein scheme 
- conv. order 2H-1/2 



mod milslein schems 
-conv. order 2H-1/2 



10 10" 10 10" 

stepsize 



Figure 4. Equation (172I) : pathwise maximum error vs. step size for four 
sample paths for H = 0.7. 



we set 

y x ■= C, ■■= r. 

Local considerations. Consider a time < T < T and for any y G Q^([0, Tq\), define z = 
T To (y) as the unique process in Q«([0,T o ]) such that z = y and (5z) st = Jsti^iv) dx). 
If y,y e Q*([0,T ]) with {yo,y$) = (yo,Vo) = (a,(r(a)), and if z = T To (y),z = T To (y), 
then some standard differential calculus easily leads to 

W*\ Q X MTo])] < c x {1 + T^[y; Q^O, T ])] 2 } , (73) 

and 

M\z-~z- Q x K ([0,T })} 

< c x T W[y-y; Q*([0,T ])] {l +M[y; Q%([0, T ])] 2 + Af[y; Q*([0, T ])] 2 } , (74) 

with Ce = c(l + ||x|| 7 + ||x 2 || 27 ) for some constant c > 1. Now set T = (4c 2 )~ 1//( ^ 7 ~ K ' ) and 
R To = 2c x , so that, if in addition M[y\ Q*([0, T })} < R To , then by ([73]), JV[z; Q*([0, T ])] 
< R To and, if also J\f[y; Q%([0,T o })} < R To , by flTJ), 

- 5; Q:([0,T ])] < c^fo - y; ^([0,r ])] • (4c 2 )-"/^) {l + 8c 2 } . 
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Observe that 3 — 2k/ {^f — k) < for 1/3<k<7<1/2 and so 

c x (4cl)-^ {1 + 8c 2 } = {cl-^-^ + ScT 2 ^-)} < 9 (jj < 1. 

As a result, IY is a strict contraction of the following closed subset of Q^([0,T ]): 
B Z<r(a)),R To ={y£ Q%(^To}); (yo,y x Q ) = (a,a(a)) , Af[y; Q x K ([0,T })} < R To } . 

Let us denote by y T ° the fixed point of the restriction of IY to B^° a ^ Rt . 

Extending the solution. One can use the same arguments as in the previous step for the 
set 

B 2T ° 

= {ye Q:([To,2T ]); (y To ,y x Q ) = (y%,a(y%)) , A%; Q X K ([T ,2T ])} < R To } , (75) 

and this provides us with an extension of the solution on [To,2To], denoted by y 2T °. 
Repeat the procedure until [0, T] is covered, and then define 

i=l i=l 

where Nt is the smallest integer such that Nt ■ T > T. 

It is not hard to see that y is a solution to the system ffl5|) . Moreover, 

< sup A/y To ; Ql([(k - 1)T , kT })} + {1 + ||x|| 7 } ^ U[y kT ^ Q*(p - 1)T , kT })} 

k=l,...,N TQ ^ 

< R To + R To ■ N To • {1 + ||x|| 7 } < 2c x (1 + (T/T + 1) (1 + ||x|| 7 )) 

< 2c x (1 + (1 + 4 • T • c^-">) (1 + ||x|| 7 )) , 

which gives the estimate (12 lj) . The unicity of this solution is easy to prove due to ( 1741 . 
The details are left to the reader. 

6.2. Continuity of the Ito map. We shall now prove point (2) in Theorem 12.61 For 
this, let us again introduce some notation: 

Notation: If y E Q x and y e Q x for two different driving signals x, x, define 

M[y - y; Q x > x ] = M[{y, y x ) - (y, y x )\ Q x f] ■= M[y - y; CJ] + J\f[y x - y x \ 

+ N[y*-yhcn 

Local considerations. Consider a time T > 0. From the decomposition 



% - y)st = [o-(y s ) - cr(y s )} ■ (5x) st + a(y s ) ■ S(x - x) st + [y x a'(y s ) - yV(y s )] ■ x 2 



St 



+ y>'(y s ) • [x 2 , - x 2 ,] + A st ( [a{yf - a{yf] ■ Sx + a{yf ■ S(x - x) 

+ 5 \y x a\y) - yV(y)] ■ x 2 , + *(yV(iO) • [x 2 - x 2 ] )., 
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where we have used 

(5y) st = [(id-A5)(a(y) ■ 5x + {a{y)f ■ y?)] st , 
some standard computations yield 

X[y-m Qf([o,r ])] 

< c x ^y {T*M[y - y; Q x /([0,T o ])] + \\x - x|| 7 + ||x 2 - x 2 || 27 + \a-a\} 

with 

<■,,,,, = c {1 + ||x|| 7 + ||x 2 || 27 + ||x|| 7 + ||x 2 || 27 + Af[y; Q*([0, T})} 2 + A% Ql ([0, T])] 2 } , 

for some constant c > 0. Now remember that Af[y; Q%([0,T])} < Pt(||x|| 7 , ||x 2 || 27 ), as 
well as Af[y; Q%([0, T])] < Pt(\\x\\^, ||x 2 || 27 ), for a certain polynomial function P T , so 
that c Xi x, y ,y < c XtX , where c X)X > stands for a polynomial expression of ||x|| 7 , ||x 2 || 27 
and ||x|| 7 , ||x 2 || 27 . Set T = (2c X:X )~ 1 ^ K and in this way 

A% - y; Q x /([0, T ])] < 2c^ {||x - z|| 7 + ||x 2 - x 2 || 27 + \a-a\}. 

Extending the inequality. With the same arguments as in the above step, we get, for any 
k > 1, 

M[y-y; Q x /([kT , (k + 1)T„])] 

< {||a; - x|| 7 + ||x 2 - x 2 || 27 + |y fcTo - y fcTo |} 

< 2c XtX | ||x - x|| 7 + ||x 2 - x 2 || 27 + \a-a\ + T oJ2^ly ~ & QT&T , (I + 1)T ])] J 
and clS 8l result 

A% - j/; Qf([kT , (k + 1)T ])] < 2c x , x ■ e k {\\x - x|| 7 + ||x 2 - x 2 || 27 + \a-a\} 

using the discrete version of Gronwall's Lemma. 
Inequality (1221) is then a direct consequence of 

n Tq -i 

M[y-y;CU[Q,T})}< E AA[y - ([A;To, (A; + 1)T ])], 

fc=0 

where iVr is the smallest integer such that Nt ■ T > T, so that Nt < 1 + T/T < 
1 + T • (2c x , £ ) K . 
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